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Abstract 

In this paper, we discuss the prompt J/ip production at e + e~ colliders via color- 
singlet and color-octet production mechanisms. The color-singlet production 
processes include 1) e + e~ — > J/ipgg; 2) e + e~ — > J/ipcc; 3) e + e~ — > qqggJ/ip 
and e + e~ — > qqgXc followed by Xc J/^l- The color-octet production pro- 
cesses include 1) e + e~ — > J/ipg; 2) e + e~ — > J/ipqq. Of all these production 
channels, we find that the color-octet contributions dominate over the color- 
singlet contributions at any energy scales. At low energies (yfs < 20GeV), the 
dominant channel is e + e~ — ► J/ipg whereas at high energies e + e~ — > J/ipqq 
will take the leading part. We also find that the energy spectrum for the 
color-octet J/ip production in process e + e~ — > J/ipqq is very soft, and the 
mean energy of the produced J/ip is only about lOGeV ~ 20GeV even at 
very high energies (e.g. at lOOOGeV). The extraction of color-octet matrix 
elements from J /if) production in e + e~ collisions is also discussed. 

PACS numberfs): 13.60.Le, 13.20.Gd 
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I. INTRODUCTION 



Since the first charmonium state J/ip was discovered in 1974, the studies of heavy quarko- 
nium states (include charmonium and bottomonium) have played an important role in ele- 
mentary particle physics. Charmonium and bottomonium are the simplest quark- ant iquark 
composite particles, which can be described by the gauge theory Quantum Chromodynamics 
(QCD). The investigation of their properties, such as the production and decays may help us 
to reveal the properties of QCD in both perturbative and nonperturbative sectors. In par- 
ticular, the Jfij) production is of special significance because it has extremely clean signature 
through its leptonic decay modes. Therefore if the heavy quarkonium production mechanism 
is clarified, J/ip triggers can be used as a powerful tool in studying other interesting physics. 

The studies of heavy quarkonium production in recent years are mainly stimulated by 
large discrepancies between the color-singlet model (CSM) predictions and the recent experi- 
mental data of CDF at Fermilab Tevatron []J . There is orders of magnitude disagreement 
between them. In the color-singlet model description for the quarkonium states, the quark- 
antiquark pairs are created in colorless configurations. In the past few years, a rigorous 
framework for treating quarkonium production and decays has been advocated by Bodwin, 
Braaten and Lepage in the context of nonrelativistic quantum chromodynamics (NRQCD) 
0. In this approach, the production process is factorized into short and long distance parts, 
while the latter is associated with the nonperturbative matrix elements of four-fermion oper- 
ators. Another outstanding feature of NRQCD is that it treats the quarkonium not simply as 
a quark-antiquark pair in color-singlet but rather a superposition of Fock states. Although 
generally the color-singlet takes a more important role in quarkonium production and de- 
cays, the other high Fock states may be dominant in some cases. Under this framework, 
one can calculate the inclusive production and decay rates to any order in strong coupling 
constants a s as well as v 2 , the relative velocity of heavy quarks inside the bound state. Just 
with this mechanism, the authors of Ref.[3-5] have successfully explained the J/i/j(ip ; ) and 
T production surplus problems discovered by CDF group. 

J ftp production in e + e~ annihilation process has been investigated by several authors 
within color-singlet model[6-ll]. Recently, Braaten and Chen have noted that a clean signa- 
ture of color-octet mechanism may be observed in the angular distribution of J / ip production 
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near the endpoint region at e + e collider []12"|. In this paper, we make a thorough discussion 



about the J/ip production at e + e~ colliders, including color-singlet and color-octet contribu- 
tions. The rest of the paper is organized as follows. In Sec. II, we discuss the J/i[) production 
color-singlet processes, including 1) e + e~ — > J/ipgg; 2) e + e~ — > J/ipcc; 3) e + e _ — > qqggJ/ip 
and e + e~ — > qqgXc followed by Xc J l^l- We carefully study the characters of each chan- 
nel. In Sec. Ill, we study the color-octet production mechanism, including 1) e + e~ — > J/ipg\ 
2) e + e~ — > J/ipqq. Comparisons of the energy scaling properties as well as energy spectrum 
of these processes are made. In Sec. IV, we give a short discussion about the determina- 
tion of color-octet matrix elements, and find that they can be extracted from e + e~ collision 
experiments at different energy regions. Finally, we briefly discuss our results and make a 
conclusion in Sec.V. 

II. COLOR-SINGLET CONTRIBUTIONS 



The leading order color-singlet contributions to direct J/ip production include the fol- 
lowing processes 

e + e~ — > 7* — > J/ipcc, (1) 
e+e- -> 7 * -> J/ip99, (2) 
e + e~ -»• 7* -> qqg* with g* -> J/ipgg, 

(3) 

The relevant Feynman digrams are shown in Fig.l. Fig. 1(a) is a quark process, Figl.(b) is a 
gluon process, Fig. 1(c) and Fig. 1(d) are gluon jets processes. All these three processes have 
been calculated separately in literatures [8- 11]. Here, we make a comparison of the relative 
weight of them function of cm. energy y/s. 



For the quark process, the calculation is straightforward, and we get JTT 



4 

da{e+e- -> J/^ee) = 32ageg < Oj ( 3 S 1 ) > fi/^ 

o^dzdx x ' 243 m 3 (1 - x 1 ) 2 {2 - z) 2 (z + x x - l) 4 ' U 

where 
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2p ■ k 2pi ■ k m 2 

z = , Xi = , r = , \p) 

s s s 

k, p and pi are the momenta of virtual photon 7*, J/ip, and outgoing parton (c quark 
in Fig. 1(a), gluon in Fig. 1(b)), respectively; m is the mass of J/ip, which equals to 2m c 
at nonrelativistic approximation; < Of ( 3 Sx) > is the color-singlet nonperturbative matrix 
element; The functions Fi are defined as 

Fi = 2(xi + z- l) 2 (si - l)[4x? - 6x 2 (z 2 - 2z + 2) 

-2xi(4;z 3 - Ibz 2 + 12z - 2) - llz 4 + 23z 3 - 24z 2 + 8z + 4], 
F 2 = [x x + z - l)[16xf + 4x 4 (5z - 4) + 4x?(3z 2 - 4z - 24) + x 2 (5z 3 + 16z 2 - 180z + 256) 

+xx(-7z 4 + 30z 3 - \32z 2 + 360z - 256) + \\z A - 55z 3 + 136z 2 - 200z + 96], 
F 3 = 8(4^ - 14x 4 + 31xf - 50x 2 + 43x x - 14) + 4z(14x 4 - 62x? + 149a; 2 - 184zi + 83) 

+2z 2 (21x'l - 130x 2 + 269xi - 176) + z 3 (29x 2 - 154x x + 177) + 4z 4 (4x 1 - 11) + 5z 5 , 
F 4 = 3{4xl + Ax\z - 4x1 ~ ~ z 3 + Az 2 - 8z + 4)(2xi + z - 2). (6) 

After integrating Eq.(Q) over Xi, and considering the fragmentation approximation (i.e., 
\fs 3> m c , r«l), the differential cross section reads as 

^(e + e _ -> J/Vcc) = 2cx(e + e~ -> cc) x D c ^j/^(z), (7) 

where D c ^j/,^(z) denotes the quark fragmentation function. Its explicit form may be found 



in Ref. 13 



As for the gluon process Fig. 1(b), from Ref. || we readily have 



da( e +e- -> J/jjgg) _ 64e 2 c a 2 <Of{ 3 Sx) > 
a^dzdxx 81 m 3 

where 

. (2 + x 2 )x 2 (2 + Xx)xx 

f(z, xx;r) = — r^y- ^ + 



r /(z,xi;r), 



+ 



(2-z) 2 (l-xi-r) 2 (2 - z) 2 (l - x 2 - r) 2 

(z-r) 2 -l 1 ■ 



1 / 6(l + r-z) 2 

— z) 2 ^ (1 — x 2 — r) 2 (l — xx — 



(1 — x 2 — r) 2 (l — X! — r) 2 (2 — z) 2 ^(1 — x 2 — r) 2 (l — Xi — r) 2 

2(l-*)(l-r) | 1- 

1 — x 2 — r)(l — Xi — r)r r> 



+ „ 2(1 "l (1 " r) , +-), 0) 



the variables z, Xi, r are defined as Eq.(|) and x-i = 2 — z — x\. 

J ftp production from gluon jets processes as shown in Fig. 1(c) and Fig.l(d) have been 
calculated in Ref. [[[(]], giving 

s 

a ( e + e - ^ qj}g *. g * _» j/^x) = J dfi 2 a(e + e- - qqg*(^))P(g* - J/^X), (10) 

in- 

where \x = m(g*) is the virtuality of the gluon. P(g* — > J/ipX) is the decay distribution 
function of virtual gluon to J/ip, which includes the contributions directly from g* — > J/ipgg 
and from El transitions of Xc (g* —* gXc followed by Xc —> J/ipj)- We can express it as 

P(g* - J/^X) = P s (g* - J/^gg) + Br( Xc - J/^l)P{g* - Xc9), (H) 

where 



lQqg < gffiS^j > r 
2437T m 3 



P s (<7* -> J/^) = ^ J / / dan/^mjr). (12) 



2^ »- 

Here the function f(z,X\] r) is defined as Eq.([|) with r = m 2 /fi 2 . The integration limits of 
X\ are 

1 



x± = -(2-z± Vz 2 -4r). (13) 

Because the g* —>■ Xcjg processes have the infrared divergence involved, which are associated 
with the soft gluon in the final states, we induce an infrared cutoff to avoid the singularities 
in P(g* — > Xcjg)- Strictly speaking, the divergences can be cancelled in the framework of 
NRQCD (see Ref. 0). Here, we follow the way of Ref. |D| by imposing a lower cutoff A on 



the energy of the outgoing gluon in the quarkonium rest frame. Then, the decay distribution 
functions for Xcj can be written as 

P to . ^) = i ^WIMi - 30' (14) 

97T n z m° 1 — r r m c 



8a 2 \R' P (0)\ 2 6r(l + r) 1 A 

P{g -> Xci£ = -7; — — i 0( 1), (15) 

97r yU^m 1 — r r m c 

8a? |ffp(0)| 2 2r(l = +r + 6r 2 ) fl ,l A 

^£ Xc2^ = ~ — : 0( !)• (16) 

97T 1 — r r m c 



As discussed in Ref. [|15|], the cutoff A can be set to m c in order to avoid the large logarithms 
in the divergent terms. 
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With the above formulas, we can evaluate the prompt J/ip production rates of color- 
singlet processes in e + e~ annihilation at any energy regions. The results are displayed in 
Fig. 2 and Fig. 3, the input parameters used in the numerical calculations are (16| (17| 



m-u, = md = m s = 0, m c = 1.5GeV, m b = 4.9GeV, a s {2m c ) = 0.26, (17) 
< Of ( 3 Si) >= 0.73GW 3 , \R' P {0)\ 2 = 0.125GeV 5 . (18) 



The branching ratios of Xcj —> J/ipl taking in the calculations are [18 



Br(xdo - 




= 6.6 x 10 






= 27.3%, 


Br(xc2 - 




= 13.5%. 



(19) 

The angular distribution and energy distribution of color-singlet J/ip production at 



CLEO have been discussed in Ref. JTTJ] . Here, we make a comparison of the relative impor- 
tance of these three color-singlet processes as a function of cm. energy yfs. The result is 
displayed in Fig. 2. The dotted line demonstrates the contribution from quark process (1), 
the dash line is from gluon process (2), the dotted-dash line is from gluon jets processes (3), 
and the solid line is the total color-singlet cross section. From this figure, we can see that 
at low energies (< 25 GeV) the gluon process dominates the other two processes, at some- 
what high energies the quark process will dominate, and at high enough energies, the gluon 
jets processes are dominant. To see the sensitivity of quark fragmentation approximation 
to interaction energy y^s, we also plot the line corresponding to the quark fragmentation 
approximation Eq.(7) in the same figure. Obviously, the quark process at high energies 
can be represented by the quark fragmentation approximation. The diagrams show that 
at -y/i > 70 GeV the difference between the complete calculation and the fragmentation 
approximation is less than 5%. Another striking result is that the process (3) is negligible 
at low energies, but at high enough energies (a/s > 200GeV) its contribution dominates the 
other two processes and grows with the energy increases. 

In Fig. 3, We display the energy distributions of quark process, gluon process, and their 
sum at y's to be 10.6 GeV, 25GeV, 50GeV^ and 100 GeV respectively. At low energies, such 
as at CLEO (-\/i = 10.6GeV^), the energy spectrum of these two processes are both flat 
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(see Fig. 3 (a)), however, at high energies the patterns of the energy spectrum of these two 
processes are distinct. The quark process is hard and the gluon process is soft (see Fig. 3(c)- 
(d)). After the relative importance of these two processes is changed with cm. energy, the 
energy distribution feature of the total cross section of these two processes is also changed 
with cm. energy. Therefore, at high enough energies, the distribution is mainly from quark 
process, and the total spectrum appears hard (see Fig.3(a)-(d)). 

III. COLOR-OCTET CONTRIBUTIONS 



The leading order color-octet contributions to direct J/ip production in e + e collision 
include the following two processes 

i) e + e~ ^7* ^qq + cc[8, 3 S x \, (20) 

ii) e + e-^7*^S + cc[8, 2S+1 Lj], (21) 

as shown in Fig.4. Here 2S+1 Lj denotes the states 1 5 , and 3 Pj, q represents the u, d, s, c 
and b quarks. 

As the first process shown in Fig. 4(a), using the factorization formalism described in Ref. 
0, we can write the differential cross section as 

da{e + e- -> qqJ/vp) = da(e + e- -> qqcc[8, 3 S x ]) < Of^S^ >, (22) 

where da represents the short distance coefficient of the process, which can be calculated 
perturbatively. < Og ( 3 Si) > corresponds to the long distance nonperturbative matrix ele- 
ment. It can be treated as free parameter or evaluated by fitting the theory to experimental 
data. The result is 

3 

da{e+e - _> qqjfr) = e y g < q^ Si ) > fi/^ 

a^ t dsids2 12 s 2 (si — m^)(s — S\ — s 2 — m 2 — m 2 ,) ' 

where 

s 1 = (k-p 1 ) 2 , s 2 = (k-p) 2 . 

In the above, k is the momentum of virtual photon 7*, p and p\ are the momenta of outgoing 
J I ip and quark q, m q is the mass of the quark and e q is its charge. The integration limits of 
S\ and S2 are 
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m 2 + m 2 q - — [s 2 (s 2 - s + m 2 ) ± A2 (s 2 , s, m 2 )\^ (s 2 , m 2 m 2 )], (24) 
zs 2 

s 2 =Am% 4 = (v^-m) 2 . (25) 



A is defined as 



The functions G,- are 



\(x, y,z) — x + y + z — 2xy — 2yz — 2xz. 



Go = 3m 6 m 2 + 4m 2 m® + 7m 4 m 4 + 2m® + s(m 6 + 4m® + m 6 + llm 4 m 2 + 16m 2 m 4 ) 
+s 2 (2m 4 + \\m 2 m 2 + 7m 4 ) + s 3 (m 2 + 3m 2 ) + Si(— m 6 — 6m 4 m 2 — 5m 4 s 



— 12m 2 m 4 — 16m 2 m 2 s — 5m 2 s 2 — 8m® — 12m 4 s — 6m 2 s 2 — s 3 ) + s 2 (3m 4 



+ 12m 2 m 2 + 8m 2 s + 12m 4 + 12m 2 s + 3s 2 ) + 4s?(-m 2 - 2m 2 - s) + 2s 4 



q - 3 ) 1 "°D 

G*i = —5m m 2 — 8m 2 m 4 — 2s (m 4 + 6m 2 m 2 + 4m 4 ) — s 2 (2m 2 + 5m 2 ) 

+sx[m 4 + 4m 2 m 2 + 4m 4 + 4s(m 2 + m 2 ) + s 2 ] + 4s 2 (-m 2 - 2m 2 - s) + 4s 1; 
G 2 = 3m 2 m 2 + m 2 s — m 2 si + 3m 4 + 3m 2 s — 2m 2 si — ss\ + 3s 2 , 

G 3 = Sl - m 2 . (26) 

A check can be performed by considering the high energy limit in this process. At 
high enough energies, the outgoing quark mass m q can be neglected. Setting m q = 0, and 
integrating over s\ and s 2 , we will obtain 



3 



a(e + e~ -> Jj^qq) e 2 q a 2 (2m c ) < Ot{ 3 S 1 ) > 



07^ 96 m 3 



{5(1 — r 2 ) — 2r lnr 



+ [2Li 2 ( T ^-)-2Li 2 ( I l-) 

- 2 ln(l + r) lnr + 3 lnr + In 2 r](l + r) 2 }, (27) 

where Li 2 (x) = — J dt ln(l — £)/t is the Spence function. The result here is completely 


consistent with that in Ref . |19| , in which the charmonium production in Z° decays through 
the similar process is calculated. 



The second process as shown in Fig. 4(b) has been calculated in Ref. |12[], giving 



„( e + e - ^ J/^g) =C ,< 0*CS a ) > +C P < 0*( 3 P a ) >, (28) 

C. = ( 29) 
6 s z m 

_ 256yr 2 aV r(l - 3r) 2 6(1 + r) 2(1 + 3r + 6r 2 ) 
9s 2 m 3 1 — r 1 — r 1 — r 



where we have used the approximate heavy quark spin symmetry relations: 

< O+fPj) >^ (2 J + 1) < OtfPo) > • (31) 

Up to now, the color-octet matrix elements < C$( 3 Si) >, < Of^So) > and < C^( 3 P ) > 
are determined only by fitting to the experimental data. < Og( 3 Si) > is extracted from 
hadroproduction process at the Tevatron[2][3] |L6||. The results are consistent with the 
theoretical anticipation of NRQCD. As for < 0$ ( 1 5'o) > and < Of ( 3 Po) >, they have been 
obtained from both hadroproduction and photoproduction, but the results from different 
processes are incompatible (see Ref. p0[), which we will give a detailed discussion in Sec. IV. 
Here, we tentatively choose the color-octet matrix elements obtained in Ref. M [IB|, which 
are consistent with the velocity scaling rules, 

< Oi^SJ > = 1.5 x 10- 2 GeV 3 , (32) 

< OfCS ) > = 10- 2 GeV 3 , (33) 

< °« (3 f o) > = 10- 2 GeV 3 . (34) 

Using the input parameters as in Sec. II, we can calculate the intermediate color-octet 
contributions to J/ip production in e + e~ collisions. ^From Eq.(23)-(26) and Eq.(28)-(34), 
the contributions of these two processes to J/ip production rates as a function of cm. energy 
\fs may be obtained as shown in Fig. 5. The dotted line is from process (ii), the dotted-dash 
line is from process (i), the solid line is the sum of them, and the dash line is the total 
cross section of color-singlet processes contributions. From this diagram, we can see that 
at low energies the dominant process is channel (ii), at higher energies (y/s > 20GeV) the 
channel (i) dominates. It is interesting to note that the color-octet contributions dominate 
the color-singlet contributions at any energy values of ^/s. So if the color-octet production 
mechanism is correct and the color-octet matrix elements are not far small with values used 
here, J/ip production in e + e~ collision mainly comes from color-octet contributions. 

Another important character of the color-octet J/ip production in e + e~ collision is its 



energy spectrum. The energy spectrum of process (ii) has been discussed in Ref. [12]. In 
this process, ip 's energy mainly fixes at the endpoint region which equals half of the cm 
energy, its width is about {^f x 500) MeV. In contrast, the energy spectrum of process 
(i) is very soft. As shown in Fig.6, the energy distribution of J/ijj production in process 
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(i) is mainly from low energy region even at high cm. energies ( > 500GeV). Numerical 
result shows, the mean energy value of J/tp from process (i) is about lOGeV ~ 20GeV at 
any energy value of yfs less than lOOOGel^. 

IV. COLOR-OCTET MATRIX ELEMENTS 



Unlike color-singlet matrix elements which associate with the quarkonium radial wave 
function at the origin in the nonrelativistic limit, and can be calculated by potential model 
T7| , the color-octet matrix elements are unknown. They can be extracted from experi- 



mental data or from lattice QCD calculations. Before lattice QCD giving out the results, 
the color-octet matrix elements are determined only by fitting the theoretical prediction 
to experimental data. As done in the literatures [2-4] [T(J [EIJ pT| , the color-octet matrix 
elements are extracted from experimental data of J/ip production in hadron collisions and 
e~p collisions. In these processes, the production mechanism is associated with the struc- 
ture function of hadrons, therefore there are still left a lot of uncertainties, e.g. the higher 
twist effects, which we can not handle clearly now. Therefore, the extraction of octet matrix 
elements from these different processes certainly exists large theoretical uncertainties, even 
their results are not consistent with each other [^0] [21]. In contrast, the mechanism of 
J/ip production in e + e~ annihilation process is much clearer than those hadron processes 
discussed above. The parton structure is simpler, and there is no higher twist effects to 
be considered, so the theoretical uncertainty is much smaller. Following, we discuss the 
possibility of extracting elements < Of ( 3 Si) >, < Of ( 1 5o) > and < Og( l So) > from e + e~ 
annihilation experiments. 

^From the results of Sec. Ill, the color-octet contributions dominate the color-singlet 
contributions at any energy regions. The total cross section of J/ip production in e + e~ 
collision is sensitive to the color-octet matrix elements. So, we can precisely extract them 
from fitting the theoretical prediction to the experimental data. Furthermore, there is an 
another important feature of color-octet J/ip production in e + e~ process, that the relative 
importance of these octet processes are distinct at different energy regions. At high energies, 
the element < Og( 3 Si) > is important, while at low energies, < Of( l S ) > and < Of( 3 P ) > 
are important. So, we can extract them separately from different energy experiments. 
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< Of ( 1 5'o) > was determined previously from hadroproduction at the Fermilab Teva- 
tron. In Ref. |16[], the authors use the gluon fragmentation approximation at high P± to fit 



the experimental data. They obtain 

< OfCSi) >= 1.5 x 10~ 2 GeV 3 . (35) 

In Ref. 0] , the authors extend out the gluon fragmentation region, and take a global fitting 
to the experimental data (at P± > 5GeV) including < Of ( 3 Si ) >, < Of( l S ) > and 

< Of ( 3 -Po) > contributions, get 

< OfCS^ >= 6.6 x 10~ 3 GeV 3 . (36) 

At LEP II energy region (s/s ~ 160GeV), the color-octet production cross section erg is 
about 6 times larger than the color-singlet process cross section <7i, and the former mainly 
comes from the contribution of < Of( 3 Si) > (the contributions from < Of ( 1 5'o) > and 

< Of ( 3 -Po) > are small enough and may be neglected, see Fig. 5). So, < Of ( 3 <Si) > can be 
extracted by precisely measuring the total cross section of J/ip production in e + e~ annihi- 
lation at this energy region. 

As for the elements < Of ( 1 5o) > and < Of ( 3 -Po) >, the situation is more complicated. In 
previous studies, the values of them were extracted from the experimental data of hadropro- 
duction at high P±, photoproduction at forward direction and fixed-target hadroproduction. 
In Ref. 0], a global fitting to all P± region data shows that at low P± boundary the theo- 
retical prediction is dominated by the contributions from < Of ( 1 5'o) > and < Of ( 3 Po) >, 
and the fitted result is 



< OfCS ) >+^< Of( 3 P ) >= 6.6 x 10- 2 . (37) 

mi 



However, the studies of photoproduction at e p collisions show that the matrix element 



values in above equation may be overestimated [20], and the authors obtain an another 
linear combination of these two elements 



< OfCS ) >+^< Of( 3 P ) >= 2.0 x 10- 2 . (38) 
mi 



Eqs.(|3?|) and ( |38| ) are incompatible. Furthermore, the fixed-target result [21J gives the same 
argument for the matrix elements values in Eq.([J7]), and gives 
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< OfCSo) > +^ < Ot( 3 P ) >= 3.0 x 10- 2 . (39) 

This problem may be further clarified in J/ip production in e + e~ collision experiment. 
Because at low energies the J/ip production dominantly comes from the color-octet 1 Sq 
and 3 Pj subprocesses (see Fig.5), and where the associated color-octet matrix elements are 
< Og ( 1 5'o) > and < Of ( 3 -Po) >• We can extract an another linear combination of these 
two elements by fitting to the experimental data at this energy regions. Furthermore, the 
coefficients in front of these two elements in the combination are different at different cm. 
energy, because the relative importance of these two subprocesses is changed with energy. So 
we can extract different combinations from different energy experiments. Here we choose two 
typical cm. energy values, i/i = 4.6GeV (at BEPC) and 10.6GeV (at CLEO), to see what 
combination of the two elements can be extracted from the experiment. At = A.6GeV, 
from Eq . (|28|) we get 

A 8 = 0.065 < OfCSo) > + < ^ (3 9 Po)> , (40) 
and at a/s = 10.6GeV, the combination is 



A 8 = 0.26 < OtCSo) > + ^ 8v : u ^ . (41) 



< Ot( 3 P ) > 



m 2 c 



^From these two equations, the individual value of these two elements may be extracted. 



V. CONCLUSIONS 



In this paper, we have calculated the direct J / ip production in e + e~ annihilation including 
color-singlet and color-octet contributions. We have studied the energy scaling properties 
as well as energy spectrum of all the production processes. The numerical result shows that 
the color-octet contributions dominate the color-singlet contributions at any cm. energy 
scales and the energy spectrum of these two color-octet processes are distinct. That enables 
us to further carefully study the properties of color-octet J/tp production in e + e~ collision 
experiment such as the angular distributions et ai. 

In this paper, we have concentrated on the J/ip production only through the virtual 
photon in e + e~ collisions. The contributions from the Z° boson at high energies should be 
included and will be considered elsewhere. 
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Because in e + e processes J/ip production has much smaller theoretical uncertainty 
THAN IN HADRONIC J /if) PRODUCTION PROCESSES, it can be used to extract the 
color-octet matrix elements precisely. At high energies color-octet 3 Si subprocess is dominant 
and the element < Of ( 3 Si) > can be extracted from experimental data in this energy region. 
At low energies color-octet l So and 3 Pj subprocesses will dominate, one can extract two sets 
of linear combinations of the elements < Og ( 1 5'o) > and < Of ( 3 -Po) > from different energy 
experiments, e.g. at BEPC and CLEO energy regions. ^From these two combination 
equations, the individual values of the matrix elements can be obtained separately. This 
really provides a strong motivation in experiment to extract the color-octet matrix elements 
in e + e~ colliders at now reaching energy. In conclusion, in e + e~ annihilation experiment, 
J/ijj production signature provides an another criterion in testing the color-octet signals and 
the NRQCD scaling rules. 
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Figure Captions 



Fig.l. Feynman diagrams for direct J /if) color-singlet production processes in e + e~ annihi- 
lation, (a) quark process e + e~ — > J/ipcc; (b) gluon process e + e _ — > J/ipgg; (c) gluon jets 
process e + e~ — ► ggt?* with gr* — > J/tpgg; (d) x c production from gluon jets in e + e~ — > ggg* 
with g* -> Xcfi 1 - 

Fig. 2. Color-singlet cross section vs cm. energy. Dotted line illustrates the quark process, 
dash line shows the gluon process, dotted-dash line comes from gluon jets, and the sum of 
the three processes is plotted as solid line. The quark fragmentation approximation is also 
shown as short-dashed curve. 

Fig. 3. Color-singlet energy distribution da/dz as function of z at different cm. energy. The 
distributions from quark process (dotted line), gluon process (dashed line) along with the 
sum of them (solid line) at (a) yfi = lO.GGeV; (b) yfe = 25GeV; (c) yfi = 50GeV; (d) 
y/s = lOOGeV. 

Fig. 4. Feynman diagrams for direct J /if) color-octet production processes in e + e~ annihila- 
tion, (a) e + e~ — > qqJ/ip; (b) e + e~ — > J/if)g. 

Fig. 5. Color-octet cross section vs cm. energy. Process (i) contribution is represented as 
the dotted-dash line, process (ii) as the dotted line, and the sum as the solid line. The 
color-singlet contribution is drawn as dashed line. 

Fig. 6. Energy spectrum of produced J /if) in color-octet process (i) at different cm. energies. 
From up to down, the curves represent ^ = lOOOGeV, 500Ge\/, 200GeK, lOOGeK, and 
50Ge\/ respectively. 
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